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Abstract. Experimental results on partial photo-ionisation cross sections of 
helium are analysed in the light of recent advances in the semiclassical theory 
of two-electron atoms. Byun et al [T] predict that the total photo-ionisation cross 
section below the double-ionisation threshold can, scmiclassically, be described 
in terms of contributions associated with classical orbits starting and ending in 
the triple collision. The necessary modifications of the semiclassical theory for 
partial cross sections is developed here. It is argued that partial cross sections 
are also dominated by the triple collision dynamics. The expected semiclassical 
contributions can be identified in the Fourier transformation of the experimental 
data. This clearly demonstrates for the first time the validity of the basic 
assumptions made in [l]. Our findings explain furthermore in a natural way 
the self-similar structures observed in cross section signals for different channel 
numbers. 
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The rich resonance spectrum of two-electron atoms below the double ionisation 
threshold has been explored both experimentally and numerically up to principle 
quantum numbers TV w 13 — 17 of the remaining one-electron atom after ionisation 
[HISlllllSllSlEllHlini HO] • Progress towards even higher values thus moving closer 
to the three-particle break-up threshold E — Q are hampered experimentally by the 
limited photon energy resolution and numerically by the high-dimensionality of the 
system. The complexity of the classical dynamics and the large number of degrees of 
freedom of the system have so far also restricted semiclassical calculations of individual 
resonances to subsets of the full spectrum and again small N values [TTl [T^] . Until 
recently it was thought that similar restrictions also apply to a semiclassical treatment 
of total and partial photo-ionisation cross sections for E < 0. The resonance density 
increases dramatically for energies approaching the double-ionisation threshold and 
individual resonances overlap and interfere leading to a strongly fluctuating cross 
section signal which decreases in amplitude towards the threshold [H [H [3] . The strong 
interaction between resonance poles can be regarded as a signature of the underlying 
chaotic classical dynamics which is also reflected in the resonance spacing distribution 
which shows a gradual transition towards that of random matrix theory approaching 
the double ionisation threshold from below [21 [6l \W\ . 

We will consider in the following partial photo-ionisation cross sections of two- 
electron atoms in the asymptotic regime E ^ 0_. Near the threshold, electron- 
electron correlation effects dominate which can be observed directly in scaling laws 
such as Wannier's celebrated threshold law for double ionisation [l3l [14] or in slow 
electron ionisation experiments taken across the threshold [15]; a cusp-like structure 
in this "zero-kinetic energy" ionisation cross section is observed which has been 
interpreted in terms of classical escape along the Wannier ridge leading again to a 
threshold law with Wannier's exponent both below and above E = [16]. Zero- kinetic 
energy spectroscopy below the double ionisation threshold does, however, not resolve 
the complex resonance structure of the three-particle compound and the experimental 
signal contains thus little information about the mostly chaotic scattering dynamics 
of the underlying classical three-body Coulomb system. 

Semiclassically, total photo-ionisation cross sections can be described in terms of 
closed orbit theory (COT). The theory was developed for systems such as hydrogen in 
external fields [iTl [18] or in the context of quantum defect theory for many electron 
atoms [m [20l [21] for which the dynamics near the origin is regular. Recently, the 
necessary modifications for a closed orbit treatment of the total photo-ionisation cross 
section for two-electron atoms have been presented in [Tj. Starting point of a COT is 
the total cross section in dipole approximation written in the form 

a{E) = -At: ahLo^{D(j3,\G{E)\D(j)^) (1) 

where (pi is the wave function of the initial bound state, D = 7r-r is the dipole operator, 
TT is the polarisation of the incoming photon with angular frequency lu, and G{E) is 
the Green function of the system at energy E = Ei + huj; furthermore a = e^ /he is the 
fine-structure constant. In the semiclassical limit, the support of the wave function 
shrinks to zero relative to the size of the system reducing the integration in ([T]) to 
an evaluation of the Green function at the origin in the limit i? — s- 0_. Writing the 
Green function in semiclassical approximation |17[ I18[ I22| thus leads to a summation 
over contributions from classical trajectories starting and ending at the origin. 

For two-electron atoms, the origin r = (ri,r2) = represents the point where 
both electrons reach the nucleus simultaneously, that is, all three particles collide. Note 
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that we work in the infinite nucleus mass approximation, that is, the position of the 
nucleus is fixed at the origin. The presence of such triple collisions demands a careful 
re-evaluation of COT in the light of the three-body dynamics near the origin. The 
triple-collision itself forms a non-regularisable singularity of the classical equations of 
motion, that is, trajectories ending in a triple collisions can not be continued through 
the singularity. Trajectories coming close to a triple collision become extremely 
sensitive to initial conditions and nearby orbits approaching the collision point can be 
scattered into arbitrarily large angles. The triple collision singularity is in that sense 
infinitely unstable. This is in contrast to binary collisions which can be regularised by 
a suitable space and time transformation such as described in [23| leading to a smooth 
phase space flow in the vicinity of the collision. 

A semiclassical treatment of photo-ionisation starting from ([T]) needs to take 
into account classical trajectories beginning and ending near the three-body collision 
R = where i? = (r^ -I- r^)-'^/^ is the hyper-radius; the set of closed orbits usually 
employed in COT, namely those emerging out of and returning exactly to the triple 
collision point i? = 0, are infinitely unstable and give a vanishing contribution to 
semiclassical expressions. These closed triple collision orbits (CTCO) act, however, 
as guiding centres for phase space regions leaving and returning to the triple-collision 
region. Taking the semiclassical limit 0_ is equivalent to Rq ^ in appropriately 
rescaled coordinates where Rq characterises the size of the initial wave function 0^. 
By considering these limits carefully, it has been shown in [1] that the amplitude of 
the fluctuations in the total photo-ionisation cross section decays with a power-law 
according to cr/; oc \E\^ with predicted exponent 



with Z, the charge of the nucleus. One obtains, for example, /i — 1.30589... for 
helium with Z = 2. The exponent can be obtained as a combination of stability 
exponents of the triple-collision singularity and differs from Wannier's exponent which 
describes double ionisation processes. The Fourier components of the fluctuations can 
furthermore be associated directly with CTCO's. This behaviour has been confirmed 
by quantum calculations in coUinear helium, the restricted three-body Coulomb 
problem where the dynamics takes place along a common axis [1] . 

In this letter, we analyse the experimental data on partial photo-ionisation cross 
sections for helium presented in [H [5] (and shown in Fig. [T^), in the light of the 
theoretical predictions [l] . The experiment was conducted using synchrotron radiation 
with a resolution of 4 meV revealing partial cross sections up to energies 78.85 eV above 
the ground state and reaching ionisation channels of the order iV w 13. The kinetic 
energy of the outgoing electron was measured which yields information on the state of 
the remaining He'*' - ion; for experimental details, see [4]. Ref. [5] notes in particular, 
that the cross sections crjv for different channels show similar patterns up to an overall 
shift in energy, see for example the data for TV = 4 and A'' = 5 in the region E < —0.02 
au in Fig. [T^; (in Fig. [1] and throughout the paper, atomic units (au) are employed). 
The phenomenon could be reproduced in R-matrix calculations published in the same 
article. In Fig. [TJd, the fluctuating part of the experimental partial cross sections are 
shown after subtracting numerically a smooth background contribution. 

In the following, we will argue semiclassically that CTCOs also dominate the 
fluctuations in partial ionisation cross-sections and show that collision orbits can 
actually be detected in the experimental data. This leads naturally to a semiclassical 




(2) 



Evidence of triple collision dynamics in partial photo-ionisation cross sections of helium 4 




* "J/ !V.;> V y ■ " V -V-^-V-- ■vV^..•~■•^.;'^A..■^■-•■-A^V.■•.«,/V■"■•.• N = 7 

: '"■^:v^;^i^^^-:'vv■"v••■- n = 9 



N 


= 4 


N 


= 5 


N 


= 6 


N 


= 7 


N 


= 8 


N 


= 9 


N 


= 10 



-0.03 



-0.02 



-0.01 







r< = i2 

N-=-ll 

N=10 

N = 9- 



r 






N = 4- 
_i 



N: 

■N=:6 

x:5 



■N;=8- 
= 7 



-0.03 



-0.02 



-0.01 







E(au) 



Figure 1. a) Partial cross sections as measured in [4] [5]; (the jump in 
the signal for N = at E —0.02 au is due to experimental reasons; 
the signal for N = A was in addition shifted down towards N = 5). b) 
The fluctuating part of the cross sections for different channel numbers 
obtained from a) after subtracting a smooth background contribution. 



explanation for the similarities seen in the cross-sections ajy for different N. The 
experimentally observed slight shifts in energy between the patterns in different partial 
cross sections can then be explained in terms of phase differences between direct and 
indirect contributions as discussed in more detail later. The experimental data clearly 
show a decrease in the amplitude of the fluctuations as i? — *■ 0_; the exponents 
characterising the mean power-law decay of the fluctuations can at present not be 
extracted from the data with sufficient accuracy to allow a comparison with theoretical 
predictions; they will thus not be considered here. A detailed theory of threshold laws 
for partial cross sections will be presented elsewhere |24| . 

We start by expressing the partial cross sections ct^t in terms of the retarded Green 
function G{E), from which a semiclassical theory can be developed. The outgoing 
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solution of the inhomogcncous Schrodinger equation 

(E ~ H)\x >^ D\^, > (3) 

is written as 

\x>^GiE)\^,> with 1$, >=i?|0, >, (4) 

and using the notation introduced after Eq. ([T]). In the asymptotic hniit R oo 
which corresponds for fixed E' < to either ri or r2 — > oo, the Schrodinger equation 
(PI is homogeneous and separable. One obtains asymptotically for the inhomogeneous 
solution \x > in the ionisation channel (A^A) 

< r^NXlx >=< ri7VA|G(S)|$, >cx /+(ri) for n ^ cx) (5) 

where A = {LMI1I2} denotes the quantum numbers describing angular momenta and 
N is the principal quantum number of the remaining ion; furthermore, f^{ri) is the 
outgoing Coulomb function 

Z — 1 

cxp[i(fcArri + — Inri)] (6) 

Kn 

with kjsi = \/2{E — En), the momentum of the outgoing electron. After writing the 
cross section as flux through the surface R = const, taking the limit R ^ 00, and using 
the asymptotic form of the inhomogeneous solution (j4|) as given in ([5|) , we express the 
partial cross section crjv in the form [25j 

(TAT = Vctata = V27rawfcAr lim \G{riNX;E)f (7) 

Z J Z J ^QQ 

A A 

where G{riNX; E) =< riNX\G{E)\^i > and cttv becomes independent of ri. 

The matrix elements G{riNX; E) of the Green function can in semiclassical 
approximation be described in terms of classical trajectories starting at R < Rq near 
the origin and reaching ri ^ 00 with fixed energy En of the hydrogen like atom 
and fixed quantum numbers A, thus determining the kinetic energy cn ~ of 
the escaping electron. For details see j24| ; a similar treatment applied to transport in 
quantum wires with fixed channel numbers can be found in |26j . We will distinguish 
direct escape leading from the initial region R < Rq directly to ionisation and indirect 
contributions from trajectories entering into a predominately chaotic phase space 
region before ionisation, that is, 

G(riiVA; E) = Ga^MNX; E) + G,nd{riNX; E). (8) 
The cross section unx can now be written in terms of a smooth background part 
and a fluctuating contribution according to (Jnx = Cjva + "^tva 'v^here 

al,^ =2TTaujkN\Gd;r{riNX-E)\'' (9) 
= A^aujkN 3? [Gl„{riNX- E)G,„drr{riNX; E)] . (10) 

In (|10p . the term |Gi„dirP is neglected as it contains contributions from pairs of 
indirect trajectories which give rise to lower order corrections in the asymptotic limit 
\E/En\ « 1. 

We turn now to a brief discussion of the classical dynamics in two-electron atoms. 
By introducing the scaling transformation [23| 

r = f/|£;|; p^VWlP; S = S/^/\El L = L//^, (11) 
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one studies the dynamics at fixed energy E ~ —1 where f , p denote the scaled 
coordinates and momentum and L is the total scaled angular momentum. The region 
R < Ro containing the initial state 0; shrinks according to R^ — l^^l-Ro ~^ for 
E' ^ 0_, that is, the initial conditions approach the triple collision. Likewise, the 
kinetic energy of the outgoing electron diverges according to cat = l\E\ for fixed N . 
Expressing the angular momentum in scaled coordinates, we have L ^ as £' — > 
for fixed L; the classical dynamics takes place closer and closer to the zero angular 
momentum manifold. It can therefore asymptotically be characterised by trajectories 
in the invariant subpace L = which has only three degrees of freedom [T^ . The part 
of the dynamics contributing to the semiclassical Green function G = Gdir + Gmd 
is in scaled coordinates formed by trajectories starting closer and closer to the triple 
collision i? = as (x i?o — ^ 0. The triple collision singularity i? = itself has a 
non-trivial structure and is equivalent to the phase space dynamics a,t E = Q. We will 
here describe only those features of the dynamics near R = which are important 
for understanding the semiclassical approximations; for more details see [T^ \Tl\ 
Most initial conditions near the triple collision will give rise to trajectories leading to 
immediate ionisation of one electron. Direct orbits escaping with kinetic energy eN 
will contribute to the direct escape term Gdir', an example of such a trajectory (dashed 
line) is shown in Fig. Only a fraction of the phase space near R = can enter a 
chaotic scattering region; these orbits all move out along the so-called Wannier orbit 
(WO), the trajectory of coUinear and symmetric electron dynamics with ri = r2. 
An example of such a indirect orbit is shown in Fig. [2K (thick full line). Likewise, 
trajectories inside the chaotic phase space can only leave this region after approaching 
the triple collision along the WO. Escaping from the triple-collision region towards 
ionisation with asymptotic kinetic energy ~ ejv/|-E| ^ oo as E' ^ 0_ can be 
achieved only by coming closer and closer to the triple collision. The point of closest 
approach, i?7v, of trajectories escaping with kinetic energy e^v vanishes for fixed N like 
Rpf (X \E\. The indirect part of the Green function Gind{N) is thus semiclassically 
described in terms of orbits starting and ending (in scaled coordinates) closer and 
closer to the triple collision before escaping with kinetic energy cm- These trajectories 
approach a proper CTCO asymptotically. The phase space regions contributing to 
the partial Green function Gindir{N) can thus be characterised in terms of CTCO's 
in the same way as the total cross section [1] . Triple collision orbits only occur in the 
so-called eZe space [12], a collinear subspace of the full three body dynamics where 
the two electrons are on opposite sides of the nucleus [23]. As i?o, — ^ 0, only orbits 
coming close to the eZe space can start and return to the triple collision and they will 
do so in the vicinity of a CTCO. The dynamics in the eZe space is relatively simple as 
it is conjectured to be fully chaotic with a complete binary symbolic dynamics. The 
symbolic coding of a trajectory is here defined as: 

— : electron 1 collides with the nucleus - ri = 0; 
+ : electron 2 collides with the nucleus - ^2 = 0; 
c : triple - collision. 

The completeness of the symbolic dynamics implies that there is exactly one CTCO 
for every finite binary symbols string; the shortest is the WO with code cc. In Fig. 
[2b, some short CTCOs are shown; note that orbits whose symbol code is related by 
the operation -|- — are mapped onto each other by the particle exchange symmetry 
Ti ^ r2 and are thus equivalent. 

CTCOs have been shown to be important for the total cross section [I] giving 
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Figure 2. Triple collision dynamics and closed triple collision 
orbits: a) trajectories starting near the triple collision point 
ri = r2 = and leading to ionisation with fixed energy of the 
escaping electron; dashed (green) line: direct path; thick full (red) 
line: indirect path, here close to the CTCO cH — c; b) some short 
CTCOs together with binary symbol code. 



the main contributions in a modified semiclassical closed orbit treatment of photo- 
ionisation. The fact that CTCOs also enter partial cross section here is due to the 
above-mentioned escape mechanism which is special to two-electron atoms. Swarms of 
indirect trajectories starting at R ~ Rq, following a similar path in the chaotic region 
and then escaping with fixed kinetic energy e^v can for ^ be linked to a single 
CTCO together with a direct-escape orbit. Likewise, the phase space region leading 
from R < Rq to direct escape can be represented by a single trajectory. We thus write 
Gd^r ~ Ao(iVA)e''^«(''i^^'^'-''^«('-i^^'^)"/2 . (^2) 

r , ~ JSo{riN\;E)-ii^oiriN\;E)7T/2 

X ^ Aj{NX;E)e'^^'-^^-"'^''/\ (13) 
CTCO 

Here 5*0, vq are the action and Maslov index along the direct orbit, respectively, and 
Sj, Vj are those along a CTCO starting from and ending in the triple collision. g| The 
initial wave function <f>i and the dynamics in phase space leading to direct orbits is not 
sensitive to the threshold energy E = Q. The pre-factor Aq will thus pick up a smooth 
energy dependence and can be treated as constant at the threshold. The singular 
behaviour of the classical dynamics near the triple collision for CTCO contributions 
leads to a universal scaling behaviour of the terms Aj as £^ — > 0_ as described in [T]; 
the semiclassical analysis suggests the same scaling law as for the total cross section, 
that is, 

Aj{N\■E)^aJ{NX)\E\^' (14) 

for ^ 0_ and fixed N with exponent given by Eq. (g]). From d?]), dTU]), (O, (O and 
p4p. we obtain an approximation to the fluctuating part of the partial cross sections 

X We note that extra phase contributions arise due to the stationary phase approximation leading 
to the fixed kinetic energy condition in the outgoing channeh these phases are not energy dependent 
and are here absorbed in the pre-factors. 
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where the pre-factor Aj{N) are of the form Aj{N) = 

z ~ note, that the Aj{N^ are in general complex valued and will add 

extra (energy independent) phases. 

The semiclassical treatment above suggests that a Fourier transformation of a)^ 
in terms of the variable z would reveal peaks at the actions of classical CTCOs. Such 
an analysis has been carried out using the experimental data from [4l [5] setting the 
origin of the energy scale at the double ionisation threshold, that is, 79.014 eV above 
the ground state of helium. After subtracting a smooth background by fitting a low- 
order polynomial to the experimental data, the fluctuating signal as shown in Fig. 
[TJ) is transformed using the Lomb - algorithm [55], a Fourier transform technique 
particularly useful for unevenly spaced data sets. (Note that the experimental data 
arc taken at uniform steps in E which leads to an increasing step-size in z.) The peak 
heights in the experimental partial cross sections are roughly of the same size for fixed 
energy independent of A^. The corresponding Fourier signals thus show structures 
of similar size, see inset in Fig. [3] to enhance the resolution, we averaged over the 
transformed cross sections for iV = 4 — 8 for which data sets in energy intervals of the 
same length are available. The resulting averaged Fourier signal is shown in Fig. [3] 
together with the values of the actions of short CTCOs such as depicted in Fig. [SJd. 

Despite the fact that the energy range available is rather limited, there are clear 
correlations between the positions of the pronounced peaks and the actions of CTCOs; 
(note, that the high-resolution Fourier signal presented in [T] for the total cross section 
was obtained using a much larger energy range - up to = 50 - in the coUinear eZe 
model). The Fourier signal shown in Fig. [3] is a clear indication that photo-ionisation 
cross sections of three dimensional helium are indeed dominated by the triple collision 
dynamics in the coUinear eZe subspace. The similarities in the partial cross sections 
noted in [5] , arc now naturally explained in terms of the semiclassical expansion ; 
the dominant terms contribute asymptotically with the same phase - the actions of 
CTCOs - which will produce similar modulations in the overall signal for different A^. 
The slight energy shifts observed in the cross section pattern when comparing signals 
for nearby A^ values as can be seen Fig. [1] are caused by energy independent phases in 
the Aj (N) 's entering Eq. ([T5|) ; a detailed analysis will be presented in [21] . 

Due to the degeneracy of the c + c and c — c orbits related by particle exchange 
symmetry, the peak associated with the code -|- is about twice as high as the WO - 
peak. A similar phenomena was observed for the total cross section in [T]. Otherwise, 
the semiclassical amplitudes decrease in general exponentially with the length of the 
symbol code which can also be seen in the Fourier data. Furthermore, a peak at an 
action smaller than that of the shortest CTCO, the WO with action S/2tt = 3.5, can 
be observed in Fig. [3] Wc note, that the peak is at S'/27r = V2, which is the action 
of the asymptotic threshold orbit where the outer electron escapes with zero kinetic 
energy leaving the inner electron at a total energy En = — 1. The peak may thus stem 
from non-perfect cancellations between actions for direct and indirect contributions 
especially near the channel thresholds. It is a feature special to partial cross sections 
and is not expected to be observed in total photo-ionisation signals. 



in the form 




(15) 
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Figure 3. The Fourier spectrum of the fluctuating part of tlie 
partial cross section in Fig. [1] averaged over A'' = 4 — 8; the actions 
of CTCOs are marked on the axis together with their respective 
binary symbol code. Inset: individual Fourier signals. 



In conclusion, wc show that the fluctuations in the partial photo-ionisation cross 
section in helium below the double ionisation threshold are semiclassically dominated 
by contributions from closed triple collision orbits. This explains naturally the 
similarities observed in different partial cross sections and is a clear evidence that 
both partial and total cross sections are dominated by the low-dimensional collinear 
eZe dynamics [T] . It is at present not possible to extract threshold laws such as stated 
in Eq. ((M]) from the data; a detailed account of the theory including the iV-dependence 
of the amplitude terms due to the direct and indirect contributions will be presented 
elsewhere 
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